NCR GAL LIBRARY FUNCTIONS 


The GAL programs included in this listing are "library" functions 
which may be used in conjunction with the GAPSYS software package of 
the GAPP PC Development System (NCR45GDS1-MS). The intention is to 
aid the user in developing GAL ‘(GAPP Algorithm Language) programs for 
the GAPP processor (NCR45CG72). These programs serve two purposes: 
they provide the user with some examples of how GAL subroutine 
programs may be written and they provide a collection Of library. 
functions which may be used for applications programming. Note that 
these programs have been written as "general purpose" routines. In 
many cases, it is possible to optimize (modify) the code for a 
specific application to obtain superior performance over these 
standard functions. 


The following functions have been implemented as GAL subroutines: 


absval Find absolute value of 2's complement number 

acctc Accumulate 2's complement numbers 

addsm Add signed magnitude numbers 

addtc Add 2's complement numbers 

addum Add unsigned magnitude numbers © 

compare Compare unsigned or 2's complement numbers 

dilate8 Dilate binary image in eight directions 

dividesm Divide signed magnitude numbers 

dividetc Divide 2's complement numbers 

divideum Divide unsigned magnitude numbers 

edge Extract edge of gray scale image 

erodes Erode binary image in eight directions 

fswitch Conditional reversal of two inputs (two output 
multiplexor ) 

meadian3 3 by 3 median filter 

median31 3 by 3 median filter, line scan version 

multsm. | ‘tude numbers 

multtc © numbers 





multum Multiply unsigned magnitude numbers 


mux Conditional selection of two inputs 
(multiplexor) using simple algorithm 

mux2 Optimized multiplexor using "tricky" algorithm 

sort Add data to a sorted list 

shift Shift in any of 8 directions 

skel Extract topological information from an image 


of an object and produce a skeleton that 
represents the original topology 


sm2tc Signed magnitude to 2's complement converter 
subsm Subtract signed magnitude numbers 

subtc Subtract 2's complement numbers 

tc2sm 2's complement to signed magnitude converter 
thresh Flag all pixels above a given threshold 

xsobel Find gradient in x direction using sobel xform 
ysobel Find gradient in y direction using sobel xform 


A more detailed description of a function may be found in the header 
section of that function's source listing. The header section will 
give a description of the function, the parameters, additional GAPP 
memory required within the subroutine, execution times (in GAPP 
cycles), and any special notes about characteristics or usage of the 
function. 


Each of these functions is found in a file listing which has the same 
name and a file extension of ".gal". For example, the "“acctc" 
function is in a file named "“acctc.gal"”. 


SPECIAL NOTE: 


All of the signed magnitude (sm) library functions allow (and 
produce) both positive and negative zero. If this is not desirable 
for a specific application, the results of these functions may be 
checked for the negative zero condition at the end of the sm 
function. If the result is negative zero, the sign bit should be 


inverted. 


/* ABSVAL.GAL */ 


/* Function: 
absval takes the absolute value of the input image (a 2's complement 
format number). This function returns the an usigned magnitude number 
in the original image. 
Parameters: HE 
image A = on input: two's complement format number 
on output: the absolute value of the input 
in unsigned magnitude format 
Additional Memory Used: 
Execution Time: 
4-bit image 17 cycles 
8-bit image 33 cycles 
16-bit image 65 cycles 
Notes: 
ny 
absval (A) 
image A; /* Input image */ 
int i, 
n=size(A); 
A:n-l c:= ram ns := Q; /* get sign bit */ 
for (i = O; i < mn; itt) /* invert if sign bit = 1 */ 
{ 


A:i ew := ram; 
A:i ram := sm; /* sm = c xor ew */ 


} 
/* Add sign bit to complete 2's complement */ 


for (i = O; i < nz; i++) 
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continues on next page 
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/* ACCTC.GAL */ 


/* Function: 
= acctc is a GAL program which adds two two's complement numbers and 
* places the result into the second operand (i.e. accumulates). 
* © 
* Parameters: 
*% image A = Ist operand 
* image B = 2nd operand ~ 
* 
* Additional Memory Used: 
* 
* ‘Execution Time: 
* 4-bit images - 12 cycles 
* 8-bit images 24 cycles 
* 16-bit images 48 cycles 
* 
* Notes: 
* 1) The second operand (image B) is overwritten with the result. 
* 2) The result is truncated, leaving the least significant bits (up 
* to the size of image B) in the second operand. 
ef 
acctc(A,B) 
image A; 
image 8B; 
{ 

int bit; /* bit counter */ 

int nA = size A, /* size of 1st operand */ 

ns = size B; /* size of 2nd operand */ 

A:0 ns:=ram c:=0; 

B:0 ew: =ram; 

B:0 ram:=sm c:=cy; 

for (bit = 1; (bit < nA) && (bit < nB); bit++) 

{ 


A:bit ns:= ram; 
B:bit ew:= ram; 
B:bit ram:= sm C:= cy; 


} 
if( nA < nB ) 
{ 
for (bit = nA; (bit < nB); bit++) 
{ 
A:nA-1l ns:= ram; 
B:bit ew:= ram; 
B:bit ram:= sm Ci= Cy; 
} 
} 


/*®* pg 1 of ADDSM.GAL */ 


+ 


Function: 
addsm is a GAL program which adds two signed magnitude numbers. 


# 


Parameters: 
image A = I1st operand 
image B = 2nd operand 
image R = *Result field 


Additional Memory Used: 
This routine uses two (2) additional GAPP RAM locations. 


Execution Time: 


4-bit images 62 cycles 

8-bit images 122 cycles 

16-bit images 242 cycles 
Notes: 


1) The result area, R, must be 1 bit larger cna the 
larger of the two operands. A general form of the 
declaration could be: 


R: (size(A) > size(B) ? size(A)+1l : size(B)+1); 


+e ee + + HF HF A OH FF HF HF FH HF OF HF KF HEN 
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addsm(A,B,R) 


image A,B,R; 


{ 
image xA:1, /* used to "zero pad" list operand */ 
xB:1; /* used to “zero pad" 2nd operand */ 
int i, /* standard counter */ 
m = size(A), /* size of ist operand */ 
n = size(B), /* size of 2nd operand */ 
r= (m>n?min)+1; /* size of result field */ 


/* take two's complement of A if sign of A is 1 */ 


Az:m-1 c:=ram ew:=0; 
for(i = O; i < m-1; itt) /* invert every bit except sign bit */ 
¢ | 


A:i ns:=z=ram; 
A:i ram:=sm; 


} 
for(i = O; i < m; i++) /* add sign bit to complete two's comp */ 


A:1l ns:=ram; 
A:i ram:=sm c:=cy; 


/* pg 2 of ADDSM.GAL */ 


/* take two's complement of B if sign of B is 1 */ 


¢ 


B:n-l c:sram ew:=0; 
for(i = O07; i < n=l; i++) /* invert every bit except sign bit */ 


B:i ns:=ram; | 1 
B:i ram:=sm; 

} 

= = 0; i <n; itt) /* add sign bit to complete two's comp */ 
B:i ns:=ram; 
B:i ram:=sm c:=cy; 

} 


/* now add numbers using two's complement arithmetic */ 


A:0 ns:=ram c:=0; 
B:0 ew: =ram; 
R:O ram:=sm c:=cy; 


for(i = 1; i < (m<n?min); i++) 
{ 
A 
B 
R 


i ns:=ram;> 
L ew:=ram; 
1 ram:=sm c:=cy; 


} 
for(i = (m<n?min); i < rel; itt) 
{ 
1£ (m <= n) Atm-1 ns:=ram; 
else A:i ns:=ram; 
1f (n <= m) Bine-l ew:=ram; 
else B:i ew:=ram; 
R:i ram:=sm c:=cy; 


} 

/* calculate sign bit */ 
Asm-l ns:=ram; 

B:n-l ew!=ram; 

R:r-1 ram:=sm; 


/* take two's complement of R if sign of Ris 1 */ 


R:r-l ct=ram ew:=0; 
for(i = O; i < rel; itt) /* invert every bit */ 
{ : 


R:i ns:=ram; 
R:i ram:=sm; 


3 
for(i = 0; i < r-l; i++) /* add sign bit to complete two's comp */ 


i ns:=ram; 
i ram:=sm C:=cy; 


R 
R 


+ 
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addtc(A, 
image A, 


{ 
int 


} 

A:m- 
Bin- 
R: 


/* ADDTC.GAL */ 


Function: 
addtc is a GAL program which adds two two's complement numbers. 


Parameters: . 
image - A = Ist operand. 
image B = 2nd operand 
image R = _ Result field 


Additional Memory Used: 


Execution Time: 


4~bit images 15 cycles 

S-bit images 27 cycles 

16-bit images 51 cycles 
Notes: 


1) The result area, R, must be 1 bit larger than the larger of the 
two operands. A general form of the declaration could be: 


R: (size(A) > size(B) ? size(A)+1 : size(B)+1); 


B,R) 

B,R; 
7 /* standard counter */ 
m= size(A), /* size of ist operand */ 
n = size(B); /* size of 2nd operand */ 


ns:=ram c:=0; 
ewi=ram: 
ram:=sm c:=cy; 


(i = 13; i < (m<n?min); i++) 


A:i ns:=ram; 
B:i ew:=ram; 
R:i ram:=sm c:=cy; 


(i = (m<n?min); i < (m>on?min); i++) 


if (m <= n) A:z:m-l1 ns: =ram; 
else A:i ns:=ram? 

if (n <= m) Bin-l ew: =ram; 
else B:i ew:=ram; 

R:i ram:=sm cC:=cy; 


1 ns:zram; 
l ew:=ram; 
? 


(m>on?m:n) ram; =sm; 


7? 


+ 


Tg 
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/* ADDUM.GAL */ 


Function: 
addum is a GAL program which adds two unsigned magnitude numbers. 


€ 


Parameters: 
image A = 1st operand . 
image B = 2nd operand 
image R = »Result field 


Additional Memory Used: 


Execution Time: 
4-bit images 13 cycles 


8-bit images 25 cycles 
16-bit images 49 cycles 
Notes: 


1) The result area, R, must be 1 bit larger than the larger of the 
two operands. A general form of the declaration could be: 


R: (size(A) > size(B) ? size(A)+1l : size(B)+1); 


addum(A,B,R) 


image A,B,R; 


{ 


int i, /* standard counter */ 
m = size(A), /* size of lst operand */ 
n = size(B), /* size of 2nd operand */ 
r= (m>n?m:in)+1; /* size of result field */ 


A:0 ns:=ram c:=0; 
B:0 ew:=ram; 
R:0 ram:=sm ¢c:=cy; 


for(i = 1; i < (m<n?min); i++) 


{ 

A:i ns:=ram; 

Bsi ew: =ram; 

R:i ram:=sm c:=cy; 
} 


for(i = (m<«n?min):; 1 < rel; i++) 


if (m <= n) ns:=0; 

else A:i nsS:=ram; 

if (n <= m) ew:=0; 

else B:i ew:=ram; 

R:i ram:=sm C:=cy; 
} 


R:r-1 ram:=c; 


/* COMPARE.GAL */ 


/* Function: 

* compare subtracts two image variables and places the resulting carry 
* in the bit-plane flag F. If A is greater than or equal to B, then F 
* will be zero. Otherwise, F will be 1. 

* 

* To allow this subroutine to be used with the sort routine, where 

* a list of data is concatenated into one image variable, the integer 
* variables pos_A and pos_B give the offset into the A and B image 

= variables, respectively. The int variable siz gives the number of 
- bits in each entry. 

7 

* Parameters: 

* (image) A lst operand. 

* (image) B 2nd operand. 

* (int) pos A offset within list operand. 

* (int) pos_B offset within 2nd operand. 

* (image) F result flag. 

* (int) siz number of bits to compare. 

* 

* Additional Memory Used: 

* None. 

* 

* Execution: 

* 2 * siz + 2 instructions. 

* 

i A 


compare(A, pos_A, B, pos B, F, siz) 


image 


nD ok opt 


int 


A, 


? 
e 
? 


pos _A, 
pos B, 
siz; 


i; 


A:pos_A ns:=ram c:=0; 
B:pos_ B ew:=ram; 


for (i = 1; i < siz; i++) { 


A:(pos_A + i) ns:=ram c:=bw; 
B:(pos B + i) ew:=ram; 


/* DILATE8.GAL */ 


+ 


Function: 
dilate8 expands a binary image in eight directions: N, S, E, W, NE, NW, 
SE, and SW. : os 


Parameters: 
image A = ,Binary image to be dilated 
image t = temporary storage location 


Additional Memory Used: | 
This routine uses one (1) additional GAPP RAM location. 


i 


Execution Time: 
13 cycles 


Notes: 
Element designation: 
NW N NE 


W E 
SW S SE 


+ tt + eH + ee HH HH HF 


~ 


Partial products are designated by (#) 


Gdilate8 (A) 


image A; 
{ 
image t:1; 
A:0 ns:= ram ew:= ram; /* get A and begin moving */ 
/* NW and NE * / 
nS:= n ew:= w c:= l; /* Wor N > cy(1) * / 
nSsi:= ew ew:= nS c= cy; 
t:0 ram:= c ns: n ew:= e c:= il: /* NE or NW > cy(2) * / 
t:0 ns:= 0 ew:= ram c:= cy; /* (1) or (2) > bw(3) x / 
A:0 ns:= ram ew: ram c:= bw? /* get SE and SW * / 
t:0 ram:= c nst= S$ ew:i=e /* S or E > cy(4) & / 
NnS:= @w e@Wl= nS C= cy; : 
A:O Yram:= ¢ nS:= S$ eW:= W cis /* SW or SE > cy(5) * / 
A:0 ns:= 0 ew:= ram c:= cy; /*® (5) or (4) > bw(6) * / 
t:0 ew:= ram c:= bw? /* (6) or (3) > bw x / 
c:= bw; 
A:O ram:= c ns:= ram ew: ram; /* store dilated image * / 


a I A 


~ 


/* page 1 of DIVIDESM.GAL */ 


Function: | 
dividesm is a GAL program which divides two signed magnitude numbers. 


Parameters: 


image P = Dividend 

image Q = Divisor 

image R = Result field 

image S = *Scratch area 

image X = Single bit scratch area 


Additional Memory Used: 
This routine uses one (1) additional GAPP RAM location. 


Execution Time: 


4-bit images 84 cycles 

8-bit images 402 cycles 

16-bit images 1758 cycles 
Notes: 


1) The size of the result field, R, is siza(P). 

2) On exit, R contains the result (an interger) and P contains 
the remainder. . 

3) R must be cleared (contain zero) before this routine is called. 

4) S$ is the same size as the larger of P and Q. 


dividesm(Q, P, R, S) 


image P, /* dividend */ 
Q, /* divisor */ 
R, /* result area */ 
: /* scratch area */ 


image x13 /* single scratch bit */ 
int yy /* standard counters */ 
int p = size(P); /* size of dividend */ 
int q = size(Q); /* size of divisor */ 


for(j = p-2; Jj >= O; j--) 


70 c:=0 ew: =#ram; 
>j ns:=ram; | 
:0 ram:=sm c:=bw? , 


or (i = 17; 1 < q-1l; i++) 


Q:i ew:=ram; , 
1f (jti > p-2) ns:=0; 
else P:jt+i ns:ram; 
S:i ram: =sm; 

c:=bw; 


} 


/* page 2 of DIVIDESM.GAL */ 


/* invert and store borrow */ 


X: ram:=c ns:=0 ew:=c c:sl1; 
R:j ram:=sm ns:=ew: B 


for(i = 0; i+j < p-l; i++) 


{ 


fe 
:i ram:=c; 

:jtLl ewi=ram c:=0; 

>i ci=cy ew:=ram ns:=0; 
>: Cc:=bw ns:=ram; 

> J 


/* calculate sign */ 


Q:q-1 ns:=sram; 
P:p-l ew:=ram ¢c:=0; 
R:p-1l ram:=sm; 


/* page 1 of DIVIDETC.GAL */ 
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Function: 
dividetc is a GAL program which divides two two's complement numbers. 


Parameters: 
image P = Dividend 
image Q = Divisor 
image R = Result field 
image S =+Sceratch area 
image X = Single bit scratch area 
image xP = Extended sign bit for P 
image xQ = Extended sign bit for Q 


Additional Memory Used: 
This routine uses three (3) additional GAPP RAM locations. 


Execution Time; 


4-bit images 190 cycles 

8-bit images 616 cycles 

16-bit images 2188 cycles 
Notes: 


1) The size of the result field, R, is size(P)+l. : 

2) On exit, R contains the result (an interger) and P contains 
the remainder. 

3) Note that the initial values of P and Q are destroyed by this 
routine. : 

4) R must be cleared (contain zero) before this routine is called. 

5) S is the size of the larger of P+l1 and Q+l. 


+ee eee ee eH He HH HOF HH He OF OF OF OF HE OH OO OE OE OE OH OE 
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dividetc(Q, P, R, S) 


image P, /* dividend */ 
; /* divisor */ 
R, /* result area */ 
S; /* scratch area */ 
image Se ae /* single scratch bit x / 
xQ:1, /* sign bit of divisor (for sm form) */ 
xP:1; /* sign bit of dividend (for sm form) */ 
int 143 /* standard counters */° 
int p = size(P); /* size of dividend */ 
int q = size(Q); /* size of divisor */ 


/* page 2 of 


/* take two's complement of Q 


Q:q-1l c:=ram ew:=0; 


xQ:0 ram:=sc;} i* 
for(i = O07; 2 < ql: i++) {[* 
Q:i ns:=ram; : 

Q:i ram:=sm; 
} 
for(i = O07; i < q-l; i++) | hea 
{ - 
Q:i ns:=ram; 
Q:i ram:=sm c:=cy; 
} 


/*® 


/* take two's complement of P 


Q:q-l1 ram:=c; 


P:p-1 c:=ram ew:=0; 
xP:0 ram:=c; | be 


for(i = O; i < pel; i++) i 
{ 

P:i ns:=ram; 

P:;i ram:=sm;> 
} 
for(i = O07 i < pl; i++) [= 
{ 

P:i ns:=ram; 

P:;i ram:=sm C:=cy; 
} 


/* 


P:p-l ram:= c; 

/* now divide magnitudes */ 
for(j = p-1; j >= 0; j--) 
c:=0 ew: =ram; 


ns:=ram; 
ram:=sm c:=bw; 


0:0 
P:j 
$:0 
(i = 


for 1; i < q; i++) 


{ 
Q:i ewt=ram; 
1f (j+i > p-2) ns:=0; 
else P:j+i ns:=zram; 
S:i ram:=sm;> 
c:=Dw; 


DIVIDETC.GAL */ 
if sign of Q isl */ 


copy sign bit into new sign bit location * / 
invert every bit */ 


add sign bit to complete two's comp */ 


use Old sign bit location for overflow */ 


if sign of P is 1 */ 


copy sign bit into new sign bit location * / 
invert every bit */ 


add sign bit to complete two's comp */ 


use old sign bit location for overflow */ 


/* page 3 of DIVIDETC.GAL */ 


/* invert and store borrow */ 


X: ramt=c ns:20 ew:=c c:el: 
R:j ram:=sm ns: =ew; 


for(i = O; it+j < pz 1++) 
{ 
S:i ew:=ram ¢c:=0; 
c:=bw; 
S:i ram:=c; 
P:j+i ew:=sram c:=0; 
S:i c:=cy ew:=ram ns:=0; 
X: c:=zbw ns:sram; 
P:j+ti ram: =c; 


} | 
/* calculate sign */ 


xQ:0 ns:=ram; 
xP:0 ew: =ram c:2#0; 
R:p ram:=sm; 


/* take two's complement of R if sign of R is 1 */ 


Rip c:=ram ew: =0; 
for(i = O; i < p; itt) /* invert every bit except sign bit */ 
{ 

R:i ns:=ram; 

R:i ram: =sm; 


} 
for(i = O; i <= p; itt) /* add sign bit to complete two's comp */ 


R:i ns:=ram; 
R:i ram:=sm c:=cy?; 


NCR GAL Library Functions 


continues on next page — 


/* page 1 of DIVIDEUM.GAL */ 


/* Function: | 

. divideum is a GAL program which divides two unsigned magnitude numbers. 
* 

* Parameters: 

* image P = Dividend 

* image Q = Divisor 

* image R = Result field 

* image S = ,.Scratch area 

* image X = Single bit scratch area 

x 

* Additional Memory Used: 

* 

* This routine uses one (1) additional GAPP RAM location. 

* 

* Execution Time: 

* 4-bit images 138 cycles 

* 8-bit images 516. cycles 

* 16-bit images 1992 cycles 

* 

* Notes: 

* 

* 1) The size of the result field, R, is size(P). 

* 2) On exit, R contains the result (an interger) and P contains 
* the remainder. 

* 3) R must be cleared (contain zero) before this routine is called. 
* 4) S is the same size as the larger of P and Q. 

ay 


divideum(Q, P, R, S) 


image P, /* dividend */ 
Q, /* divisor */ 
R, /* result area */ 
S; /* scratch area */ 
{ 
image X:1; /* single scratch bit */ 
int 14; /* standard counters */ 
int p= size(P); /* size of dividend */ 
int q = size(Q); /* size of divisor */ 


for(j = p-l; j >= 0; j--) 
{ 


70 c:=O0 ew: =<ram; 
:j ns:=ram; 
>Q ram:=sm c:=bw; 


or (i = 13; 1 < q;z itt) 


Q:1 ew: =ram; | 
if (j+i > p-1) ns:=0; 
else P:j+i ns:=zram; 
S:i ram:=sm; 

c:=bw; 
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/* invert and store borrow */ 


X: ram:#c ns:#0 ew:=sc c:21; 
R:j ram:=sm ns:=zew; 


for(i = O; it+tj < pz i++) 
{ , s 
S:i ews=ram c:20; 
c:=bw; 
S:i ram:=c; 
P:j+i ew:=ram c:=0; 
S:i c:=cy ew:=ram ns:=0; 
X: ci=bw ns: =ram; 
P:j+i ram:;=c; 


/* 


as 


Function: 


edge 


This routine 


each 


/* page 1 of EDGE.GAL * / 


is a gal program which extracts the edge from an image. 


uses the sobel transform to find the gradient at 


pixel. The gradient magnitude is then calculated and the 


result is compared to a threshold value. Any pixel which has 

a gradient magnitude greater than the threshold value is considered 
to be part of the edge. The output image contains the original 

gray scale values at those pixels determined to be part of the 


edge, 


‘Parameters: 


image 


image 
image 
image 
image 
image 
image 


int 


Additional Memory 
This routine 


Execution Time: 
4-bit image 
8-bit image 
15-bit image 


Notes: 


the rest of the pixels have a zero gray scale value. 


A = Gray scale input image (i.e. unsigned 
magnitude number) 

X # Gradient in x-direction 

Y = Gradient in y-direction 

MIN = Minimum of X and Y 

MAX = Maximum of X and Y 

MAG = Approximate gradient magnitude 

R = Result image (all pixels contain 
zero values except those on the edge) 

T = Threshold value 


Used: 
uses one (1) additional GAPP RAM location. 


330 cycles 
534 cycles 
891 cycles 


1) The maximum input image size is 15 bits because of lack of 
available GAPP RAM. 

Larger input images may be operated on with modifications 
to this code: This routine may be broken into two or more 
routines which do not require that all of the images above 
exist at once. For example, X and Y do not need to exist at 


2) 


the same t 


ime as MAG, so X and Y could be temporary (automatic) 


images which only exist until MIN and MAX have been calculated. 


edge (A, X, Y, MIN, MAX, MAG, T, R) 


image A, 


int 


xX, 
XG 
MIN, 
MAX, 
MAG, 
R; 
T; 


/*® 


/* 
ii 


A is the gray scale input */ 

Holds X gradient */ 

Holds Y gradient */ 

Holds min(X,Y) */ 

Holds max(X,Y) */ 

Holds the approximate gradient magnitude */ 

R is the gray scale output */ 

T is trial constant for performing threshold */ 
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int i, 

n = size(A)+3; 
image F:i; 
ysobel(A, Y); 
xsobel(A, X); 
absval(Y); 
absval(X); 


3 


/* subtract the absx.- absy to get flag */ 


> Ge ns := ram c := Q; 
Ys: ew := ram; 
c := bw? 

for (i = 13; i < n;z itt) 
{ 

X:i ns := ram; 

b are ew := ram; 

c := bw; 

} 
Fs ram t= ¢; 
/* use flag to get min and max values 
mux2(X, ¥, F, MIN); ’ bol 
mux2(Y, X, F, MAX); /* 
/* MAG = MAX + MIN/2 */ 
MIN:0O c $= ram? /* 
for (i = 0; i<n- 1; itt ) 
{ 

MAX: i ns z= ram; 

MIN: itl] ew := ram; 

MAG: i ram :2 sm c :2 CY; 
} 
MAX:n-1 ns := ram ew := 0; 
MAG:n-1l ram := sm c := Cy; 
MAGin ram :2 C; 


af 


MIN = X*F + Y*®~F = min(x,Y) */ 
MAX = Y*F + X*"~F = max(X,Y) */ 


round by adding lsb */ 


/* store overflow bit */ 


/* Threshold the unnormalized Gradient Magnitude (GRAD) */ 


thresh(MAG, T, F); 


/* 


Flag (F) contains the binary edge image. 


R is produced below 


* which restores the gray-scale component to the edge image. The 
* flag is anded with the input image (A) to extract the gray scale 


ia 4 


ae 
A:0 
for (i = Q; 
{ 


ns := ram; 

ew := ram c 3= 0} 
i < size(A); i++) 
A:itl c 
R:i ram 


22 cy ew := ram; 
22 c oc 3:20; 


* values of the edge. The edge is stored in the output image (R) 
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/* ERODE8.GAL */ 


/* Function: 
. erode8 shrinks a binary image in eight directions: N, S, E, W, NE, NW, 
* SE, and SW. 
Ie . 
% Parameters:. 
* image A = Binary image to be eroded 
* image t «= -temporary storage location 
* 
* Additional Memory Used: 
* This routine uses one (1) additional GAPP RAM location. 
te 
% Execution Time: 
* 13 cycles 
* Notes: 
* 
* Element designation: 
* 
* NW N NE 
* W E 
* SW S SE 
* 
* Partial products are designated by (#) 
ate 
erodes (A) 
image A; 
{ 
image t:1; 
A:0 ns:= ram ew: ram; /* get A and begin moving a 
/* NW and NE 
ns:= n ew:= w c:= 0; /* Wor N > cy(1) - 
NSs:= ew ewWw:= nS C:= cy; 
t:0 ram:= c ns:= n ew: eae c= 0; /* NE or NW > cy(2) * / 
t:0 ns:= 0 ews= ram c:= cy; /*® (1) or (2) > ey(3) * / 
A:O0 ns:= ram ew:s ram c:* cy; /* get SE and sw x / 
t:0 ram:= c ns:= S$ ew: e c:= 0; /* S$ or E > cy(4) % / 
nS:= ew ew: nS C:= cy; 
A:O0 ram:= c ns:= S ewt= Ww c:= 0; /* SW or SE > cy(5) * / 
A:0 ns:= 0 ew:= ram c:= cy; /*®* (5) or (4) > cy(6) x / 
t:0 ew:= ram c:= cy; /* (6) or (3) > cy * / 
c:= cy; 
A:O0 ram:= c;3 /* store eroded image * / 


* 
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/ 
fswitch(A, pos_A, B, pos_B, outl, pos_1l, outz2, pos_2, F, 
image A, 
outl, 
out2, 
int pos_A, 
pos B, 
pos _l, 
pos 2, 


/* page 1 of FSWITCH.GAL */ 


Function: 


fswitch uses the flag stored in the bit-plane F to select between 


two input images. 
is copied to outz2. 
copied to outl. 
“multiplexor" algorithm: 


R= ((A* B)&F)*B 


a 


If F is zero, then A is copied to outl, and B 
Otherwise, A is copied to out2, and B is 
The method used here is a variation of the 


where R is the output, A and B are inputs, and F is the flag used to 


select between the inputs. 
implemented here are 


outl = ((A ~ B) & F) ~ A 
out2 = ((A ~ B) & F} ~*~ B 


The boolean equations which are 


To allow this subroutine to be used with the sort routine, where 
a list of data is concatenated into one image variable, the integer 
variables pos_A and pos_B give the offset into the A and B image 


variables, respectively. 


Likewise, the integer variables pos_1 and 


pos 2 give the offsets into the image variables outl and out2, 


respectively. 
entry. 

Parameters: 
(image) A lst input. 
(image) B 2nd input. 
(int) pos A offset within lst 
(int) pos B offset within 2nd 


lst output. 
2nd output. 


(image) outl 
(image) out2 


(int) pos_l offset within lst 
(int) pos 2 offset within 2nd 
(image) F flag input. 

(int) siz number of bits to 


Additional Memory Used: 
1 RAM location. 


Execution Time: 
7 * siz + 1 instructions. 


S12: 


input. 
input. 


output. 
output. 


switch. 


The int variable siz gives the number of bits in each 


siz) 


/* 
ny 
/* 
a 


int e /* 


image T:1; yas 


Load the flag into EW 


F: ew:=ram> 


Loop for each bit-plane in the input variables. 


for (i = O; i < siz; 


A:(pos_A + i) 
B:(pos B +i) c 


i ram:zsm ns:zram c:=0; 


A:(pos_A + i) 


outl:(pos_1 + 1) 


B:(pos_B + i) 
out2:(pos 2 + 1) 
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i++) { 


ns: =ram; 


>; =ram; 


c:=cy ns:=ram; 


ram: =SM; 


ns: 2#Tam; 
ram: =sm; 


standard integer counter 


used for temporary storage 


/* 
{/*® 
/* 
/* 


/* 


/* 


id 
sd 


Load 
Load 


if 
po 
ne 
if 
if 


Lt 
if 
Lif 


3 MH Mh ff fi oy 


bit i from A 
bit i from B 


oltre O*+rrFOrO 


hr 


*» 


“« 


* 


= 


s 


s 


ba] 


? 


sm = A*~B 
sm = !(A*~ 
cy = !(A™B 
cy = 0 
'(1(A°B)*B) 
sm = O”°B = 


'C'(A™B)°A) 
sm = O°“A = 


B) 
) 


it ! 


* 


wld 
al 


x 


ad 


+ 
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Function: 
median3 is a GAL program for a 3 x 3 median filter. Uses the "shift" 
subroutines to move the input data, held in the image variable "input", 
through the neighborhood PE's. As each datum arrives at a PE, the 
"sort" subroutine is used to add that datum to a sorted list of values 
held in the image variable "data". After all data movement is 
complete, the middle value is picked out of the sorted list and copied 
to the “output” image variable. 


Parameters: 
(image) input input data. 
(image) data storage for sorted list of data. 
(image) output result. 
Additional Memory Used: 
None. 
Execution: 


4 input data 1492 instructions. 
6 input data 2180 instructions. 
8 input data 2868 instructions. 
10 input data 3556 instructions. 


Notes: 
1) The image variable data must be at least as big as 
9 * size(input). 
2) The image variable output must be at least as big as 
size(input). 


median3(input, data, output) 


image input, 


/* 
oF 


/* 


ys 


/* 
my 


data, ' 
output; 


int = 


Copy the central value into the sort list. 
for (i = 0; i < size(input); i++) { 


input:1 c:=ram; 
data:i ram:=c; 


Aquire and sort the N neighbor value. 
shift s(input, input); 
sort(input, data, 1); 
Aquire and sort the NW neighbor value. 


shift _e(input, input); 
sort(input, data, 2); 


/x 
my 


/* 
P| 


/* 
*7 


/* 
a & 
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Aquire and sort the W and SW neighbor values. 


for (i = 3; i < 5; i++) { 
shift _n(input, input); 
sort(input, data, i);, 


Aquire and sort the S and SE neighbor values. 


for (i = 53; i < 7; itt) { 
shift w(input, input); 
sort(input, data, i); 


Aquire and sort the E and NE neighor values. 


for (1 = 7; i < 9; i++) { 
shift _s(input, input); 
sort(input, data, i); 


Copy the median value to the output image variable. 


for (i = 0; i < size(input); i++) { 
data:(4 * size(input) + i) c:=ram; 
output: i ram: =c; 


i i i i i i 


“~, 
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* Function: 


median3l is a line scan version of the median3 program. The input 
data is assumed to be located in three image variables: input, linel, 
and line2, which correspond to three adjacent lines of the 
two-dimensional data which‘is being filtered. "input" is the last 
line to be input, and so it is the southernmost line. "line2" is the 
northernmost line. One pass through the routine will produce one 
line of the filter output. Therefore this routine should be 
repeatedly executed (by placing the routine inside a while(1) loop 
for example). ' 


Parameters: 
(image) input southernmost line of input data. 
(image) linel middie line of input data. 
(image) line2 northernmost line of input data. 
(image) temp storage for temporary working variable. 
(image) data storage for sorted list of data. 
(image) output result. 


Additional Memory Used: 


size(input) RAM locations. 


Execution: 


4 bit input data 1490 instructions. 
6 bit input data 2178 instructions. 
8 bit input data 2866 instructions. 


Notes: 


1) The image variables input, linel, line2, and temp are assumed to 
all be the same size. 

2) The image variable data must be at least as big as 
9 * size(input). 

3) The image variable output must be at least as big as 
size(input). 


median3l(input, linel, line2, temp, data, output) 


image 


int 


/* 


input, 
linel, 
lined, 
temp, 
data, 
output; 


i; 


* Copy the central vaiue into the sort list. 


"7 


for (i = 0; i < size(input); i++) { 


lineli:i c:=ram; 
data:i ram: =c; 


/* 
af 


/* 
a 


/* 


ad 


/* 


a 


{x 
if 


/* 


*/. 


/* 


a 4 


he 
if 


/* 
ay 
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Aquire and sort the W neighbor value. 
shift_e(linel, temp); : 
sort(temp, data, 1); 

Aquire and sort the E neighbor value. 
shift_w(linel, temp); 

sort(temp, data, 2); 

Sort the N neighbor value. 


sort(line2, data, 3); 


Aquire and sort the NW neighbor value. 


shift e(line2, temp); 
sort(temp, data, 4); 


Aquire and sort the NE neighbor value. 
shift_w(line2, temp); 

sort(temp, data, 5); 

Sort the S neighbor value 

sort{input, data, 6): 


Aquire and sort the SW neighbor value. 
shift_e(input, temp); 

sort(temp, data, 7); 

Aquire and sort the SE neighbor value. 
shift_w(input, temp); 

sort(temp, data, 8); 

Copy the median value to the output image variable. 
for (i = 0; i < size(input); i++) { 


data: (4 * size(input) + i) c:=ram; 
output:i | ram: =c; 


+ *¥™ 
ba 


uf 
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Shift the bottom two lines up one position. 
for the next pass through the algorithm. 


for (i = 0; i < size(input); i++) { 


s 


/*® 
* Copy linel to lineZ. 
=] 
linel:i c:=<ram; 
line2:i ram: =c; 
/* 
* Copy input to linel. 
at f 
input:i c3=ram; 
linel:i ram: =cC; 
} 


This prepares the data 
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Function: 
multsm is a GAL program which multiplies two signed magnitude numbers. 


Parameters: ; 

' dmage A = Multiplicand 
image B = Multiplier 
image C # Result field 
image X «= ‘Scratch area 


Additional Memory Used: 
This routine uses one (1) additional GAPP RAM location. 


Execution Time: 
4-bit images 47 cycles 


8-bit images 271 cycles 
16-bit images 1295 cycles 
Notes: 


1) The size of the result field is at least m+n-1l, 
where n and m are the sizes of A and B respectively. 


/ 
multsm(A,B,C) 
image A, /* multiplicand */ 
B, /* multiplier */ 
Cc; /* result area */ 
{ 
image Riis /*®* one bit scratch area */ 
int h,i, /* standard counters */ 
n = size A, /* size of multiplicand */ 
m = size B; /* size of multplier */ 


B:0 ew:=ram; 


for (i = 0; i < n-1; i++) 
{ 
A:i ns:=ram c:=0; 
c2=cy; 
Cz:i ram:=c;3 


+ 
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for (h = 13 h < m-1; h++) 


} 


A:0 ns:=ram c:=0; | ‘ 
B:h ew:=ram; 

C:h ns:=ram ew:=0 c:=cy;> 
Cz:h ram:=sm c:=cy?; 

A; ram:=c; 


7 


for (i = 1; i < n=23 i++) 
{ 
A:i nst=ram c:=0; 
B:h ew:=ram;: 
X: @ews=ram c:=cy?; 
C:ith ns:=ram;: 
Csith ram:=sm c:=cy; 
: ram:=c> 


n-2 nsi=#ram c:=0; 
h ew:=ram; 

ewizram C:=cy; 
Csnth-2 ns!:=ram; 
Csnth-2 ram:=sm ¢c:=cy; 
C:nth-1 ram:=c; 


} 
A 
B 
X 


Az:n-1l ns:=sram; 
B:m-1l ew:=ram c:20; 
C:n+m-2 ram:=sm; 
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Function: 
multtc is a GAL program which multiplies two two's complement numbers. 


Parameters: 
image X # Multiplicand 
image Y = Multiplier 
image Z = Result field 


image zsgn = Sign of result 
image karry: = Carry bit for add 


Additional Memory Used: 
This ‘routine uses two (2) additional GAPP RAM locations. 


Execution Time: 


4-bit images 160 cycles 

8-bit images 504 cycles 

16-bit images 1768 cycles 
Notes: 


1) The size of the result field is at least nx + ny, 
where nx and ny are the sizes of X and Y respectively. 
2) X and Y are modified. 
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multtc(X, Y, Z) 


image X, /* multiplican */ 
ae /* multiplier */ 
Z; /* result */ 
{ 
int bit_x, bit_y, bit_z; /* bit counters */ 
int nx, ny, nz; /* sizes of X,Y, & Z */ 
image zsgn:1; /* sign of result(Z) */ 
image karry:1; /* storage for carry during add */ 


nx = size X; 
ny = size Y; 
nz = size Z; 


/* get sign bit of result = x sign xor y sign */ 
X3nx-1 ns:= ram; 


Yiny-l ew: ram  c:= 0; 
zsgn:0 ram:= sm; 
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/* compliment X if sign of X = 1 */ 


Xi:nx-l1 c:= ram ew:= 0; 
for( bit_x = 0; bit_x < nx; bit_x++) 
{' ¢ 


X:bit_x ns:= ram ew: 0; 
X:bit_x ram:= sm; 


} | , 
for( bit_x = 0; bit_x < nx; bit x++) 


X:bit_x ns:= ram ew:= 0; 
X:bit_x rYram:= sm- c:= cy; 


} 

/* compliment Y if sign of Y = 1 */ 

Yiny-l cz: ram ew:= 0 ; 

for( bit_y = 0; bit_y < ny; bit_y++) 

{ 
Yibit_y ns:s= ram ew:= O; 
Y:bit_y ram:= sm; 

} 


for( bit_y = 0; bit_y < ny; bit_y+t+) 
{ 


¥Y:bit_y ns:= ram ew:= 0; 
Y:bit_y ram:= sm c:= cy; 


} 
/* copy ist partial product into Z */ 


Y:0 ns:= ram; 
for(bit_x = 0; bit_x < nx; bit_x++) 
{ 


X:bit_x  ew:= ram c:= QO; 
C:= cy; 

Z:bit_x  Yram:= c; 

} 

c:= 0; 

for(bit_z = nx; bit_z < nz; bit _z++) 
Z:bit_z ram:=c; 

/* multiply X & Y( ie. shift and add) 

for(bit_y = 1; bit_y < ny; bit_y ++) 

{ 
for(bit_x = 0; bit_x < nx; 
{ 
bit_z = bit_x + bit_y; 
Y:bit_y ns:= ram; 
X:bit_x ew:= ram c:= 0; 
Z:bit_z ns:= ram c:* cy; 
karry:0 ew:= ram; 
Z:bit_z Yram:= sm c:= cy 

karry:0 rams= ¢;. 

} . 


/* get sign of X > c */ 


/* comp bit if negative > sm */ 


/*®* add */ 


/* get sign of Y > c */ 


/* comp bit if negative > sm */ 


/* add */ 
as 
bit_x++) 
/* X && Y > cy */ 
/* add */ 
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/* compliment Z if zsgn = 1 */ 


zsgn:0 c:= ram ew:= 0 ; 

for( bit_z = 0; bit_z < nz; bit _z++) 

{ 
Z:bit_z ns:= ram ew:= OQ; 
Z:bit_z ram:* sm; 


} 


for( bit_z = O; bit_z < nz; bit_z++) 

{ 
Z:bit_z ns:= ram ew:= 0; 
Z:bit z Tam:= sm c:= cy; 


/* get sign of Y Dec */ 


/* comp bit if negative > sm */ 


/* add */ 
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/* page 1 of MULTUM.GAL */ 


Function: 
multum is a GAL program which multiplies two unsigned magnitude 
numbers. 
Parameters: 
image A = Multiplicand 
image B = «Multiplier 
image C -s#= Result field 
image X = Scratch area 


Additional Memory Used: 
This routine uses one (1) additional GAPP RAM location. 


Execution Time: 


4-bit images 82 cycles 

8-bit images 354 cycles 

16-bit images 1474 cycles 
Notes: 


1) The size of the result field is at least m+n, 
where n and m are the sizes of A and B respectively. 


/ 
multum(A,B,C) 
image A, /* multiplicand */ 
; /* multiplier */ 
C; /* result area */ 
{ 
image X31; /* one bit scratch area */ 
int h,2; /* standard counters */ 
n = size A, /* size of multiplicand */ 
m= size B; /* size of multplier */ 


B:0 ew:=ram; 
for(i = O; i < n; i++) 
{ 
A:i ns:=ram c:=0; 
c:=cy; 
C:i ram:=c; 
for (h = 1; h < m; h+t+) 
ns:=ram c:=0; 


0 

h ew: ram; 

h ns:eram ew:=0 c:=cy; 
hr 
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for(i=1; i < n-1l; i++) 
{ 
A:i ns:=ram c:=0; ‘ 
B:h ew:=ram; 
X: ews=ram c:=cy; 
C:ith ns:=ram; 
C:ith ram:=sm ci=cy; 
' KX: ram:=c; 


} 
Asn-l ns:2ram c:=0:> 
B:h ew:=ram; 

X: ew:s=ram cC:=cy; 
C:nth-1 ns:=ram; 
C:n+h-1 ram:=sm ¢:«cy; 
Crnth 


Nn ram: =C; 


/* MUX.GAL */ 
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Function: 
mux is a GAL program which copies one of two input images into 
an output image depending on a flag bit. The flag determines 
whether A is copied into R or B is copied into R. The algorithm 
implemented is: , 


R:i = (Asi && F:) [| (Bsi && !F:) 
Parameters: | | 
image A = Starting address of A operand 
image B = Starting address of B operand 
image R = Starting address of the Result 
image F = Flag address 
Additional Memory Used: 
Execution Time: 
4-bit images 25 cycles 
8-bit images 49 cycles 
16-bit images 97 cycles 


Notes: 


1) size(A) = size(B) = size(C) 
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mux(A, B, F, R) 


image A, B, F, R; 


{ | 
int b A /* Standard counter */ 


F:0 ns:=ram; 
for (i = 0; i < size(A); i++) 


{ 
Asi  ewseram c:=0; | /* F && A > cy * / 
B:i ew:=ram cC:=cy; /* get B % / 
Rii ram:=c c:=0; /* 'F && B > bw * / 
R:i ns:=0 ew:=ram c:=#bw; /* () || () > bw * / 
F:0 nsi:=ram c:=bw; /* flag for next bit */ 
R:ii ram: =¢C; /* resuit > R x / 

} 


Ee Nes 


+ 


+e ee FF HF FF OF OF OF OE OE OH OO OO 


+t Ft FF H HF OF 


™ 


/* MUX2.GAL */ 


Function: 
mux2 is a GAL program which copies one of two input images into 
an output image depending on a flag bit. The flag determines 
whether A is copied into R or B is copied into R. The algorithm 
implemented is: 


R:i = ((Asi “ Bri) && F:) * Beil 


This routine is FASTER than the mux routine, but a little more 
more difficult to follow. 


Parameters: 
image A = Starting address of A operand 
image B = Starting address of B operand 
image R = Starting address of the Rasult 
image F = Flag address 


Additional Memory Used: 


Execution Time: 


4-bit images 21 cycles 

8-bit images 41 cycles 

16-bit images 81 cycles 
Notes: 


1) size(A) = size(B) = size(C) 


mux2(A, B, F, R) 


image A, B, F, R; 


{ 
int i; /* Standard counter */ 
F;0O ew:=ram; /* Load F */ 
for(i = 0; i < size(A); i++) 
{ 
A:i ns:=ram; /* Load A */ 
Bi:i c:=ram; /* Load B */ 
R:i ram:=sm ns:=ram c:=0; /* sm = (A"B) if F=0, sm = 1(A7B) if F=l */ 
B:i c:=cy ns:=ram; /* cy = 0 if F=0, cy = !(A°B) if F=] */ 
R:i ram:=sm; /* sm = (O°B) {i.e. sm = B} if FzO, 
sm = !(!(A°B)°B) {1.e. sm = A} if F=1 * / 
} 
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/* page 1 of SHIFT.GAL */ 


Function: 
shift XX is a collection of subroutines to shift an image variable 
one position in any of eight directions. 


¢ 


p © shift direction: 
n north 
s south : 
8 east 
W west 
ne northeast 
nw northwest 
se southeast 
sw southwest 
Parameters: 
(image) sre source image variable. 
(image) dst destination image variable. 


Addition Memory Used: 
None. 


Execution time: 


XX: no. of instructions: 
n, S, 8, W 2n + 1 
ne, nw, se, sw 2an + 3 


where n = size(src) 


Notes: 
1) The image variable dst must be at least as large as size(src). 


shift_n(srce, dst) 


image src, /* source image variable */ 
dst: /* destination image variable */ 
int a 
src:0 nst=ram ew: =0; /* ns = sre:0, ew unused * / 
nsi=sS; | /* ns = sre:0 shifted N * / 


for (i = 1; i < size(src); i++) { 


sre:i ns:=ram c:=ns; /* oc sre:i - 1 shifted N, 
ns = sre:i a? 4 


dst:i- 1 ram:=c ns:=s c:=0; /* ns = sre:i shifted N */ 


} 


dst:(size(src) - 1) ram: =sm; /* sm = ns | * / 
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shift_s(sre, dst) 


image Src, /* source image variable */ 
dst; /* destination image variable */ 
{ 
int 3 
srco:O0 ns:=ram ew:=0; , /* ns = srce:0, ew unused 
ns:=n; /* ns = sre:0 shifted s 
for (i = 1; i < size(src); i++) { 
sre:i ns:=ram c:=ns; /* c¢ = sre:i - 1 shifted S, 
ns * sre:i 
dst:i- 1 ram:=c ns:=n c:=0; /* ns * sre:i shifted s 
dst:(size(srce) - 1) ram: =sm; /* sm = ns 
} 
shift_e(src, dst) 
image src, /* source image variable * / 
dst; /* destination image variable */ 
{ 
int 13 
src:0 ew:=ram ns:=0; /* ew = srce:0, ns unused 
ew: =w? /* ew = srce:0 shifted E 
for (i = 1; i < size(src); i++) { 
srco:i ew:=ram Cc:=ew? /* c= sre:i ~- 1 shifted E, 
ew = srecsi 
dst:i-l1 ram:=c ew:=w c:=0; /* ew sre:i shifted E 
dst: (size(src) - 1) ram: =sm; /* sm = ew 


mi 
=] 


=} 


wis 


af 


a 
mf 


af 
df 


ad 
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shift_w(sre, dst) 


image__— src, /* source image variable */ 
dst; /* destination image variable */ 


int 13 


srco:0 ew:=ram ns:=0; 
ew:=e; 


for (i = 13 i < size(src); i++) { 
sreii ew:=ram Cc: =ew; 


dst:i - 1 ram:=c ew:=ze c:20; 


/* ews src:0, ns unused */ 
/* ns = sre:0 shifted N * / 


/* c= srcti - 1 shifted w, 
ew = sreo:i */ 


/* ew = sre:i shifted N 7 a 


dst:(size(src) - 1). ram: =sm; /* sm = ew * / 
} 
shift_ne(src, dst) 
image src, /* source image variable m/ 
dst; /* destination image variable */ 
{ 
int a 
src:0 nS:=ram; /* ns = sre:0 * / 
ns:=S; /* ns = src:0 shifted N x / 
sre:l ew:=ns ns: ram; /* ew = srce:0 shifted N, 
ns = srce:l x / 
ew: =w ns:=sS; /* ew = src:0 shifted NE, 
ns = srce:1 shifted N * / 


for (1 = 2; i < size(src); itt) { 
sre:i c:=ew @eWw:i=ns nst:=z=ram; 


dst:i - 2 ram:=c eW: sw ns:2s; 
c:=ew ews=ns ns: =0; _ JR 


dst:(size(sre) - 2) ram:=c ew:=w c:=#0; /* 
dst: (size(sre) - 1) ram:=sm; /* 


/* oc = sre:i - 2 shifted NE, 


ew = srco:i - 1 shifted N, 

ns = sre:i x / 
/* ew = sre:i - 1 shifted NE, . 

ns = src:i shifted N a 


c = MSB - 1 of sre shifted NE, 
ew = MSB of sre shifted N * / 


ew = MSB of src shifted NE */ 
sm = ew =/ 
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shift_nw(sre, dst) 


image 


int 


sre; 


ste: 


for 


dst: 
dst: 


src, /* 
dst; {* 


b 
0 ns:=ram; 
ns:=*sS; 


1 ew:=ns ns:sram;> 


ew.=e@ nSi=s; 


2; i < size(src); i++) { 
Ci=ew eW:=ns ns:=ram; 


(i = 
sre:i 


dst:i - 2 ram:=C @ew:=e ns:=s; 


C:=ew @ewr=ens ns:«Q0; /* 


(size(src) - 2) ram:=c ew:=e c:=20; /* 
(size(src) - 1) ram:ssm: /* 


shift_se(sre, dst) 


image 


{ 


int 


sre: 


Src; 


for 


ast: 
dst: 


source image variable */ 
destination image variable */ 


/* ns = sre:0 * / 
/* ns = srco:0 shifted N va ae 
/* ew = src:0 shifted N, 

ns = sre:l * / 
/* ew = src:0 shifted NW, 

ns = src:l shifted N */ 


/* co sreo:i - 2 shifted NW, 


ew = src:i ~- 1 shifted N, 

ns = sre:i * / 
/* ew = sre:i - 1 shifted NW, 

ns = sre:i shifted N * / 


c » MSB - 1 of sre shifted NW, 
ew = MSB of src shifted N * / 


ew = MSB of sre shifted NW */ 
sm «= ew * / 


ero, /* source image variable */ 
dst; /* destination image variable */ 
i; 
0 ns:=ram; /* ns = src:0 * / 
ns:=n> /* ns = sre:0 shifted s x / 
l @wi=ns ns:=ram; /* ew = sre:0 shifted Ss, 
ns = sre:1 & / 
ew:=w ns:=n; /* ew = srce:0 shifted SE, 
ns = sre:l1 shifted s */ 


(i = 2; i < size(src); i++) { 


sro:ii C:=eW @ew:2ns ns:i=ram; 
dst:i - 2 ram:=c ew:=w ns:=n; 
>=ew ewr=ns ns:=0> . /* 


(size(src) - 2) ram:=c ew:=w c:=0: | 
(size(sre) - 1) ram:=sm: tk 


/* oc src:i - 2 shifted SE, 


ew = sre:i - 1 shifted Ss, 

ns = sre:i * 
/* ew = sre:i - 1 shifted SE, 

ns = src:i shifted s x / 


Cc = MSB - 1 of sre shifted SE, 
ew = MSB of sre shifted S * / 


ew = MSB of sre shifted SE */ 
sm = ew * / 
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shift_sw(srce, dst) 


image src, /* source image variable * / 


dst; /* destination image variable */ 
{ 
int = Hi 
src:0 ns:*ram; ‘ /* ns = sre:0 x / 
ns:=n;} /* ns = sre:0 shifted S * / 
sre:l eWw:=ns ns:=ram; /* ew = sreo:0 shifted Ss, 
ns = sre:l * / 
ew: 26 ns:=n; /* ew = src:0 shifted SW, 
ns = src:1 shifted $ */ 


for (i = 2; i < size(srce); i++) { 


sreo:i Cc:sew ew:=nS ns:=ram; /* o#= sre:i - 2 shifted sw, 
ew = src:i - 1 shifted S, 
ns = sre:i * / 

dst:i- 2 ram:=c ew:se ns:=n; /* ew = src:i - 1 shifted Sw, 


ns = src:i shifted S * / 
ci=ew ew:=ns ns:=0; /* c= MSB - 1 of sre shifted Sw, 
| ew = MSB of src shifted Ss * / 


dst: (size(src) - 2) ram:=c ew:=e c:=0; /* ew = MSB of sre shifted SW */ 
dst:(size(srce) - 1) ram:«sm; /* sm = ew * / 
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continues on next page. 
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/*Skeletonization Program. This program extracts topological information from 
arbitrary images to produce a skeleton of one spot thick lines that preserve 
the topology of the original object. */ . 


/*The algorithm is implemented as follows: 

Ls A Cleaning operation which eliminates spurious points in the image 

II. A Connecting operation which fills in all connecting points in the 
image 

III. A subsequent Cleaning operation to eliminate connecting points without 
neighbors in the local 4-neighborhood 

Iv. A series of Erosion operations that terminate when the image contains 
only connecting points and endpoints. */ 


/* Global Variable Declarations */ 


image 1I10:0:5; /* I/O image * / 
image X:1; /* Threshold image * / 
image SKELETON:1; | /* Working image * / 
image END:1; /* Endpoint image * / 
image CONNECT:1; /* Connecting point image */ 


/* Main Program */ 
main () 
int i; 


image SCRATCH:1; 


IO: 5 c:= ram; 
SCRATCH:0 ram:= c; 
X:0 rams= c; 


for (i120; i<5; i++) 
{ 
TO:i+l o:=ram; 
TO:1 ram: 2c; 
} 
CLEAN IMAGE (X); 
CONNECT POINT (SKELETON); 
CLEAN IMAGE (SKELETON); 


SKELETON:O NS :2 RAM C ;:2 OQ; /* Erode image in North direction */ 
NS := S EW := NS; 


CONNECT: 0 C := CY EW := RAM NS := 0; /* Preserve connecting points */ 
C := BW; 


SKELETON:0O RAM :2 C;3 
CONNECT POINT (SKELETON); 


CLEAN IMAGE (SKELETON); 
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SKELETON:0 EW := RAM C :2 Q; 
EW := W NS := EW; 


CONNECT: 0 C := CY EW := RAM 


> =BW; 
SKELETON:O RAM := 


Ce 
CONNECT POINT (SKELETON); 
CLEAN_IMAGE (SKELETON); 
ENDPOINT (SKELETON); 


CONNECT:0 NS := RAM C := 1; 


END: 0 EW := RAM; 


SKELETON:0 C := CY EW := RAM NS 
SCRATCH: 0 RAM := SM NS := RAM; 


SCRATCH: 0 NS : =RAM; 
) C:=0; 
C:=0; 


while( GOLOW); 

{ 
ERODE NORTH (SKELETON); 
CONNECT POINT (SKELETON ); 
CLEAN_IMAGE (SKELETON); 
ENDPOINT (SKELETON ); 
ERODE EAST ( SKELETON ); 
CONNECT POINT (SKELETON ); 
CLEAN IMAGE (SKELETON); 
ENDPOINT (SKELETON }; 
ERODE SOUTH (SKELETON); 
CONNECT POINT (SKELETON ); 
CLEAN_IMAGE (SKELETON); 
ENDPOINT (SKELETON); 
ERODE WEST (SKELETON); 
CONNECT POINT (SKELETON); 
CLEAN IMAGE (SKELETON); 
ENDPOINT (SKELETON); - 


CONNECT: 0 NS := RAM C := 


END:0 EW := RAM; 


SKELETON:O C :=# CY EW :s 


SCRATCH :0 RAM := SM NS 
SKELETON:0O c:=zram; 
I0:5 ram:=c; 

C:20; 


} ‘ 
C:=0; 


/* end of MAIN */ 


e¢ 


/* Erode image in East direction */ 


NS :2= 0; /* Preserve connecting points */ 


/* The operation is complete when the 
skeleton consists of only connecting 
points and endpoints */ 


/* Execution terminates when the 
skeleton consists of only 
connecting points and endpoints */ 


:= O; 
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/Xemerseeerssrzrerereszesas sk / 


/* SUBROUTINES for SKEL.GAL*/ 


« 


/* Cleaning Operation - Eliminate stray points */ 


/* Pixel neighborhood around point X: 


D X 
C 


B 


Algorithm: X = X*(A+B+C+D) 


oid 


CLEAN IMAGE (XS) 


image XS; 
{ 


image SCRATCH:1; 


XS:0 


SCRATCH: 0 


xXS$:0 


SCRATCH: 0 


xXxS:9 


SKELETON: O 


NS := RAM C := lI; 
NS := SEW := NS; 
EW := E? 

CG t= CY> 

RAM := C; 


EW := RAM NS 
C := CY; 


:= O; 


RAM := C; 


/* 
/* 


/* 
/* 


/* 


Load image into NS register */ 
Shift image North in NS; Swap 
unshifted image into EW */ 
Shift: image West in EW */ 
Register C contains (C + D) */ 


Store in RAM */ 


Shift image South in NS */ 
Shift image East in EW */ 
Register C contains (A + B) */ 
Swap contents of C into NS */ 


Register EW contains (C + D) */ 


Register C contains 
(A +B+C#+oD) */ 


Load original image into EW */ 
Register C contains 

X*(A +B+C+D) */ 

Store "cleaned" image in RAM */ 
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/*Connecting Operation - Find connecting points in image */ 


/* Pixel neighborhood around point x: 


A B C 
D Xx F 
G H I 


¢ 


Algorithm: (A + B) + (B+ C) + (C+ F) + (F +I) + (H+12) + (G+ H) + 
(D+ G) + (A+D) > 2 if point X is a connecting point */ 
CONNECT POINT (XS) 
image XS; 
{ 
image SCRATCH1:1; 
image SCRATCH2:1; 
image T:2; 
image P1:3; 
image P2:3; 
image CXT:3; 
int i; 
int n; 
XS:0 EW := RAM NS := RAM C := OQ; /*Load image */ 
EW := E> /*Shift image West */ 
SCRATCH1:0 RAM := SM; /* SM = - (A + B) (B + 
= (G + H) (H + 
EW := W NS := S; 
SCRATCH2:0 RAM := SM EW := RAM; /* SM = (A + D) - (C + 
(D = G) - (F + 
The SM output is stored in RAM 
and simultaneously loaded into 
EW := E; /* EW = - (A + D) - 
- (D + G) - 
NS := EW EW := W C:= 0; #£«°£/* Load contents of EW into NS 
EW := W; /* EW = ~ (C + F) - 
= (F + I) - 
T:0 RAM := SM C 3:= CY; /®* SM = - [(A + D) + (C + F)] 
~ [(D + G) + (F + I)]j 
gual RAM := C C ?:= 0; /* Store carry bit */ 
n = size (T); 
i = 0; 
while (i < n) 
y bee I NS := RAM EW >= RAM; 
NS := N;3 
Plii RAM := SM C := CY; 
/* a = 
SM = - [(A + D) + (C + F) + (D+G)+ (FHI 
* / - ” 
i++; 


C) 
tr) */- 


F) 
I) 


EW */ 


af 
as 


as 


- */ 


ie 
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/* subroutine CONNECT POINT continued */ 


Plin RAM := C: /* Store final carry */ 
SCRATCH1:0 NS := RAM C := O; 
NS := N;3 ; /* (NS ) =z =< - . -_ . 
- (A + B) (B + C) 
-_ = ns x / 
EW := NS NS := S&S; 
NS := S; a /*® (NS) = ~ = - 
- (G + H) (H + I) 
ma = - * / 
T:0 RAM := SM C := CY; 
/*SM = ~ - - 
- [(A + B) + (G + H)] [(B + C) + (H + I)] 
se! 
T:1 RAM :* C C := 0; 
n = size (T); 
i = Q; 
while (i < n) 
{ 
Ts NS := RAM EW := RAM; 
EW := W; 
P2:i RAM := SM C := CY; 
/* SM = = = ag 
- ((A + B) + (B + C) + (G + H) + (H + 1)j] - 
=) 
i++; 
} 
P2:n RAM := © C := Q; 
n = size (P11); 
i = Q; 
while (i < n) 
Pl:i NS := RAM; 
P2:i EW := RAM; 
CXT: 1 RAM := SM NS := RAM C :# CY; 
/* SM = [(A + B) + (A + D) + (B + C) + (C + F) + (D + G) + (F + I) + 
(G +H) + (H + I)] */ 
i++; 
} . 
EW := C C:= 1: /* Load MSB of Sum into EW: NS current- 


ly contains the next MSB of Sum */ 
CONNECT: 0 C := CY EW:=RAM NS:=0;/* C contains a "1" if point X is a 
connecting point */ 


C := BW; 
SKELETON:0 EW: =RAM; 
* CsaCy; 
CONNECT: 0 RAM := C; /* Update connecting points */ 
SKELETON:O0 EW := RAM NS := 0; f/* Preserve connecting points */ 
C := BW; 
SKELETON:0 RAM := C; | /* Update skeleton with new connecting 


points */ 


/* page 


/* Locate endpoints in image */ 


6 of SKEL.GAL */ 


/* Pixel neighborhood around point X : 


A 
D X 
Cc 


i 


Algorithm: xX is an endpoint if x*(A +B+c¢C+#+ D) = 1 */ 


ENDPOINT (XS) 
image XS; 
{ 


image T:2; 
image T1:2; 
image P1:3; 


int i: 
int n> 


NS :*= RAM EW := RAM C := O; 


NS := S$ EW := E: 
RAM := SM C : CY; 


/* SM = (C+D) */ 


RAM :* C NS := N EW := W C :# 0: 


NS := N EW := W: 
RAM := SM C :2 CY; 
RAM := C C := QO; 
i= O; 

n = size (T); 


while (i < n) 
{ 


Ts 2 
3 gt eo 
Pi:i 
i++; 
} 
P1:0 
Pl:1 


SKELETON: 0 
END:0 
SKELETON: 0 
SKELETON:0 
/* Fill in 
Pl:2 

Pi:1 

P1:0 


SKELETON: 0 


SKELETON: 0 


EW :* C C 3:2 OQ; 
EW := RAM NS :2= O 


C := CY: 
RAM := C3 
EW := RAM NS := O; 
C := BW; 
RAM := C; 


single holes - those 


EW := RAM C := OQ; 
NS := RAM; 

C := BW NS := RAM; 
EW := C C := Q 
EW := RAM NS : 


jj ~e 


0 


C := BW; 
RAM := C; 


/* SM = (A + B) */ 


/* Pl = (A+B+C#D) */ 


/* C contains P1:0 && !P1:2 x / 


C := BW; /* C contains P1:0 && !P1:1 && !P1:2 
(l.e. C = 1 iff Pl = 001) */ 
/* Now AND result with skeleton * / 
/* Store endpoints -*/ 
/* Preserve endpoints */ 


/* Update skeleton with new endpoints */ 


elements with four four-neighbors */ 


/* C contains P1:2 && !Pl1:1 * / 


C := BW; /* C contains P1:2 && fPl:l && !P1:0 
(1.6. C = 1 iff Pl «= 100) */ 
/* Now OR result with skeleton */ 
/* Update skeleton */ 
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/* Erosion Subroutines */ 
/* Erode Pixel in North Direction */ 
/* Pixel neighborhood around point X: 


x 
A 


J 


Algorithm: xX = (A * X) */ 


ERODE NORTH (XS) 


image XS; 
{ 
image ERN:1; 
XS:0 NS := RAM C := 0; 
NS := S' EW := NS; 
<= CY; 
ERN:0 RAM := C; 


/* C = (A *® X) */ 
/* Store erosion result */ 


/* In general we want to give endpoints special treatment so as to avoid 
rapid erosion of significant lines within the image. As such, endpoints 
will only be eroded when they are oriented opposite to the direction of 
erosion. This prevents an endpoint from being eroded when the direction 
of erosion is perpendicular to the orientation of the endpoint. */ 


/* Pixel neighborhood around point X: 


A Xx B 


Algorithm: X = X * (A + B) iff point X is an endpoint */ 


/* The algorithm will vary slighly depending upon the direction of erosion. 
This particular algorithm prevents an endpoint oriented in the east or 
west direction from being eroded by an erosion in the north or south 


direction. */ 


XS:0 EW := RAM C := 1; 
EW :;= E> 
NS := EW EW := W; 
EW := W; 
END: 0 C := CY NS :* RAM EW := 0; /* C = (A + B); NS contains 
endpoints in the image */ 
ERN:0 Cc s= CY :* RAM NS := O; /* C = X * (A + B)? 
EW contains eroded image */ 
BW; /* Add back unintentionally eroded endpoints */ 
sm QO; /* Preserve connecting 


C is 
CONNECT: 0 C := BW EW := RAM NS 
Cc 


SKELETON:QO RAM := C; 


points */ 
/* Update skeleton */ 
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/* Erode Pixel in South Direction */ 


/* Pixel neighborhood around point X: 


A ‘ 
Xx 
Algorithm: X = (A * X) ¥*/ | 
ERODE SOUTH (XS) 
image XS; 
{ 
image ERS:1; 
X$:0 NS := RAM C := QO; 
NS := N’ EW :# NS; 
C ss CY; 
ERS:0 2% C; 


RAM 


/* Preserve endpoints */ 
/* Pixel neighborhood around point X: 
A . 


Xx B 


Algorithm: 


X = X * (A + B) iff point X is an endpoint */ 


/* The algorithm will vary slighly depending upon the direction of erosion. 
This particular algorithm prevents an endpoint oriented in the east or 
west direction from being eroded by an erosion in the north or south 


direction. */ 
XS:0 EW := RAM C := 1; 
EW := E; 
NS := EW EW := W; 
EW := W? 
END:0 C := CY NS :s RAM EW := OQ; /* 
ERS:0 C := CY EW :* RAM NS := Q; /* 
C := BW; 
CONNECT: O C := BW EW := RAM NS := OQ; /* 
C := BW; 
SKELETON:O RAM := C? /* 


C = (A + B)?: NS contains 
endpoints in the image */ 
C=# X * (A + B)3 

EW contains eroded image */ 


/* Add back unintentionally eroded endpoints */ 


Preserve connecting 
points */ 


Update skeleton */ 
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/* Erode Pixel in East Direction */ 
/* Pixel neighborhood around point X: 
X A ‘ 
Algorithm: xX = (A * XX) */ 


ERODE EAST (XS) 


image XS; 
{ 
image ERE:1; 
XS:0 EW :s RAM C := 0; 
:2 W NS := EW; 
C := CY; 
ERE: 0 RAM := C; 


/* Preserve endpoints */ 
/* Pixel neighborhood around point X: 


A 
xX 
B 


Algorithm: xX = X * (A + B) iff point X is an endpoint */ 


/* The algorithm will vary slighly depending upon the direction of erosion. 
This particular algorithm prevents an endpoint oriented in the north or 
south direction from being eroded by an erosion in the east or west 
direction. */ 


XxS:0 NS := RAM C := 1; 
NS := N;3 
EW := NS NS := S; 
NS := S; 
END : 0 C := CY NS := RAM EW := 0; /* C = (A + B): NS contains 
endpoints in the image */ 
ERE:0 C := CY EW := RAM NS := Q; /* C =X *® (A + B); 
EW contains eroded image */ 
C := BW: /* Add back unintentionally eroded endpoints */ 
CONNECT: 0 C := BW EW := RAM NS :# QO; /* Preserve connecting 
points */ 
C := BW; 
SKELETON:0 RAM ;: C; /* Update skeleton */ 
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/* Erode Pixel in West Direction */ 
/* Pixel neighborhood around point xX: 
A X 
Algorithm: xX = (A * X) ¥*/ 
ERODE WEST (XS) . = 


image XS; 


{ 
image ERW:1; 


XS:0 EW := RAM C := QO; 
EW := E NS := EW; 
C .2= Cre 

ERW:0 RAM := C; 


/* Preserve endpoints */ 

/* Pixel neighborhood around point xX: 
A 
x 
B 


Algorithm: xX = X * (A + B) iff point X is an endpoint */ 


/* The algorithm will vary slighly depending upon the direction of erosion. 
This particular algorithm prevents an endpoint oriented in the north or 
south direction from being eroded by an erosion in the east or west 


direction. */ 
XS:0 NS := RAM C := 1; 
NS := N; 
EW := NS NS := S; 
NS := S; 
END: 0 C := CY NS :2 RAM EW := OQ; /* C = (A + B): NS contains 
endpoints in the image */ 
ERW: 0 C := CY EW := RAM NS := OQ; /* C =X * (A + B); 
EW contains eroded image */ 
C := BW; : /* Add back unintentionally eroded endpoints */ 
CONNECT: 0 C := BW EW := RAM NS := O; /* Preserve connecting 
points */ 
C := BW; 
SKELETON:0 RAM :2 C; /* Update skeleton */ 


/* SM2TC.GAL */ 


= 


/* Function: : 
" sm2tc converts the signed magnitude input image into a two's complement 
- output image. 

* 

* Parameters: ; | 

* image A = Signed magnitude format input image 

* image B = Two's complement format output image 

* ” : 

* Additional Memory Used: 

* 

x Execution Time: 

* 4-bit images 15 cycles 

* 8-bit images 31 cycles 

* 16-bit images 63 cycles 

* 

* Notes: 

* 

*% 1) The input and the output images are the same size. 

* 

* / 


sm2tc(A, B) 


image A,B; 
{ 


int. i; /* standard counter */ 
n = size(A); /* size of images .*/ 
A:n-1 c:sram ew:=0; /* get sign bit */ 
B:n-1 ram: =c; /* copy sign bit into B */ 
for(i = O; 1 < n=l; i++) /* invert bits if sign bit = 1 */ 
{ 


A:i ns:2zram; 
Bri ram:=sm; 


} 

for(i = 0; i < n; i++) /* add sign bit to complete conversion */ 
B:i ns?=ram; 
B:i Yram:=sm C:=cy; 

} 
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* 


/* SORT.GAL */ 


Function: 
The sort routine adds the contents of the image variable input to a 
sorted list in the image variable data. It assumes that data 
has already been sorted into descending order. "data" is a collection 
of values, which are assumed to be the same size as the image 
variable “input", concatenated together. That is, for size(input) = 8, 
data:0 through data:7 holds the largest value in the PE, data:8 through 
data:15 holds the next largest value, and etc. The data in "input" is 
added to the list, and the list is re-sorted. The re-sort is 
accomplished by adding "input" to the bottom of the list percolating 
the new value up in the list until it reaches its correct position. 
The integer variable “count” is the number of items already in the 
list. 


Parameters: 
(image) input input data image variable. 
(image) data sorted list of data items. 
(int) count number of items already in the sorted list. 


Additional Memory Used: 
None. 


Execution Time: 
count * (9 * n + 3) instructions, where n = size(input). 


Notes; 
1) The value of count passed to the subroutine must be at least 1. 


sort(input, data, count) 


image input, 


int 


{ 


data; 
count; 


int i, 
pos_l, 
pos_ 2; 


image F:l; 


size(input) * (count - 1); 
size(input) * count; 


pos_l 
pos _2 


compare(data, pos_1, input, 0, F, size(input) ); : 
fswitch(data, pos_1, input, 0, data, pos_1, data, pos 2, F, size(input)); 


for (i = count - 1: i > O; i--) { 


pos 1 = size(input) * (i - 1); 
pos 2 = size(input) * 1; 


compare(data, pos_l, data, pos 2, F, size(input)); 
fswitch(data, pos_l, data, pos_2, data, pos 1, data, pos 2, F, 
size(input) );~ ~ 
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F 


Function: 
subsm is a GAL program which subtracts two signed magnitude numbers. 
The algorithm is: 


A- Be2R 
Parameters: 
image A «= I1st operand 
image B = «2nd operand 
image R = Result field 
Additional Memory Used: 


Execution Time: 


4-bit images 60 cycles 

8-bit images 120 cycles 

16-bit images 240 cycles 
Notes: 


1) R must be one bit larger than the larger of A and B. 
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~~, 


subsm(A,B,R) 


image A, B, R; 
{ 


int 1 /* standard counter */ 
nA = size(A), /* size of A operand */ 
nB = size(B); /* size of B operand */ 


/* take two's complement of A if sign of Ais 1 */ 


A:nA-1 c:=ram ew:=0; 
for(i = 0; i < nA-1; i++) /* invert every bit except sign bit */ 
{ 


A:i ns:=z=ram? 
A:i ram:=sm; 


} 
for(i = O; i < nA; i++) /* add sign bit to complete two's comp */ 
{ 


A 
A 


i ns:=ram; 
i ram:=sm c:=cy; 


} 
/* take two's complement of B if sign of Bis l */ 


B:nB-1 c:=ram ew:=0; 
for(i = O; i < nB-l; i++) /* invert every bit except sign bit */ 
{ | 


B:i ns:=ram; 
B:i ram:=sm; 


} | | 
for(i = O; i < nB; i++) /* add sign bit to complete two's comp */ 
{ 


B:i ns:=ram; 
B:i ram:=sm cC:=cy; 
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/* now subtract the two's complement integers */ 


A:0 ns:sram c:=0; /* subtract LSB w/o borrow */ 


B:0 ew:=ram; 
R:0 ram:=sm c:=bw;> 


1£(nA >= nB) | "9 
for(i = 1; i < nA; itt) 


{ 
A:i ns:=ram; 
B:(i<nB ? i : nB-1) ews=ram; /* 
R:i ram:2=sm c:zbw; 
} 
else 


for(i = 1; i < nB; i++) 


A:(i<nA 2? i : nA-1) ns:=ram; /* 
B:i ew:=ram; 
R:i ram:=sm c:=zbw; 


} 


A:nA-1 ns:=2sram; /* 
B:nB-1l ew:=ram; 
R:(nA>nB ? nA : nB) ram:=sm; 


/* take two's complement of R if sign of 


R:(nA>nB ? nA : nB) c:=ram ew: =0; 
for(i = 0; i < (nA>nB?nA-1:nB-1); i++) 
{ 

B:i ns:=ram; 

B:i ram:=sm; 


} | 
for(i = O; i < (nA>nB?nA:inB); i++) 
{ 

B:i ns:=ram; 

B:i ram:=sm c:*cy; 


sign extend B for i>nB-1l */ 


sign extend A for i>nA-1 */ 


sign extend R to allow overflow */ 


Ris 1 */ 


/* invert every bit except sign */ 


/* add sign bit to complete */ 
/* two's complement conversion */ 
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/* SUBTC.GAL */ 


/* Function: 
* subtc is a GAL program which subtracts two two's complement numbers. 
* The algorithm is: 

* 

* A-B2zR ‘ 

* 

* Parameters: 

* image A = -+lst operand 

* image B = 2nd operand 

* image R = Result field 

* 

* Additional Memory Used: 

* 

* Execution Time: 

* 4-bit images 15 cycles 

* 8-bit images 27 cycles 

* 16-bit images Jl cycles 

* 

* Notes: 

& ° 

* 1) R must one bit larger than the larger of A and B. 
*/ 

subtc(A,B,R) 


image A, B, R;> 
{ 


int A /* standard counter */ 
nA = size(A), /* size of A operand * / 
nB = size(B); /* size of B operand * / 
A:O ns:=ram c:2=0; /* subtract LSB w/o borrow */ 


B:0 ew:=ram: 
R:0 ram:=sm c:=bw; 


if(nA >= nB) 
for(i = 17; i < mA; i++) 
{ 


Az:i ns:=ram; 
B:(i<nB ? i : nB-1) ew:=ram; /* sign extend B for i>nB-1 */ 
R:i ram:=sm c:=bw; 
} 
else 
for(i = 1; i < mB; i++) 


A:(i<nA ? i : nA-1) ns:=ram: /* sign extend A for i>nA-1 */ 
B:i ew:=ram; 
R:i ram:=sm c:=bw; 


} 


A:nA-1 ns:=ram; /* sign extend R to allow for Overflow */ 
B:nB-1 ew:=ram; | 
R:(nA>nB ? nA : nB) ram:<sm: 


+ 


+e te tee eH He HEH HH eH OH OH HN 
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/* TC2SM.GAL */ 


Function: 
tc2sm converts the two's complement input image into a signed magnitude 
output image. 
Parameters: ° | 
image A = Two's complement format input image 
image B = Signed magnitude format output image 


# 


Additional Memory Used: 


Execution Time: 


4-bit images 15 cycles 

8-bit images 31 cycles 

16-bit images 63 cycles 
Notes: 


1) The output image, B, must be one bit larger than the input 
input image, A. That is, size(B) = size(A)+l. 


tc2sm(A, B) 


image A,B; 
{ 


int fi /* standard counter */ 


n = size(A); /* size of input image */ 
A:n-1 :=ram ew:=0; /* get sign bit */ 
Bin ram:=c;} /* copy sign bit into B */ 
for(i = O; i < n=-1; i++) /* invert bits if sign bit = 1 */ 
{ 


Az:i ns:=ram; 
B:i ram:=sm; 


} 
for(i = 0; i < nel; it+) /* add sign bit to complete conversion */ 
B:i ns:=ram; 
B:i Yam:=sm c:=cy; 
ae ram: =c; /* store last carry */ 


/* THRESH.GAL */ 


/* Function: 
thresh compares every pixel in an image to a specified threshold 
value and reports which pixels exceed this threshold. If a pixel 
has a value greater than than the threshold, the Flag bit is set. 


Parameters: 
image A «= Input image (i.e. unsigned 
magnitude number) | 
image F = This bit indicates whether a pixel is greater | 
than the threshold value. (F = 1 if pixel is 
greater than the threshold, F = 0 otherwise. ) 
int T = The threshold value 


Additional Memory Used: 


Execution Time: 


4-bit image 14 cycles 

8-bit image 26 cycles 

16-bit image 50 cycles 
Notes: 


a f 


thresh (A, T, F) 


image A; /* input image */ 
int T> /* threshold value */ 
image — F; /* flag */ 

int i; 


/* subtract each pixel from the threshold value */ 


c := 0; 
for (i = 0; i < size(A); i++) 


/* if bit "i" of T is 1, set c equal to 1 and move previous bw to ns */ 


1£(T & (1 << 1)) 
c := 1 ns t= Gc? 


/* if bit "i" of T is 0, set c equal to O and move previous bw to ns */ 


else 
c :2 0 ons :=2 c; 


/* load image in ew, put threshold value into ns, and put previous 
borrow back into c */ | 


A:l ew := ram ns := c ci: ns; 
Cc := bw; a 


ram := C; /* store last borrow in F */ 


hy 


/* page 1 of XSOBEL.GAL */ 


/* Function: | 
xsobel takes the sobel transform of the input image in the east - 
west direction (x direction). In other words, this function finds 
gradient of the image about a pixel in the east - west direction 
by implementing the algorithm: 


Xsobel = (F + 2I1 + L) = (D + 2G +2) 
where the 3x3 neighborhood around the pixel "H" is defined as: 
D E F 
G H I 
J K IL 
Parameters: 
image A = Gray scale input image (i.e. unsigned 


magnitude number) 

image B = Sobel transform of input image for the . 
east ~ west direction (a 2's complement 
format number). 


Additional Memory Used: 


Execution Time: 


4-bit image 60 cycles 

8-bit image 104 cycles 

16-bit image 192 cycles 
Notes: 


1) Image B is three bits larger than image A. 
iW 


xsobel (A, B) 


image A, /* input image */ 
B; /* output image (y sobel of input) */ 
int i, 


n = size(A); 
/* Compute the first partial sum */ 


* 
e 


Ons := ram c := Q; 
O ew := ram ons 3: s3 
O 


A 
A 
B ram := sm C := CY; 


for (i = 13 i < nz itt) 

{ 

A:i ns := ram; 

A:i ew := ram ns :2 Ss}? 
B:i ram := sm c¢ 3: cy; 


} 
Bin ram := sm c3:= 0: /* store overflow */ 


‘ 


/* page 2 of XSOBEL.GAL */ 


/* Compute second partial sum * / 
for (i = 0; i < n+l; i++) 
{ ‘a 


B:ins := ram; 
B:l ew := ram ns :z n; 
B:i ram := sm C 3 Cy; 
} « 
Bin+l ram := sm ci := O> 
/* Compute difference */ 
for (1 = 0; i <n + 23 i++) 
{ 


B:i ew := ram: 
ew := w> 
B:i ew := ram 
ew := e; 
B:il ram := sm c 3:23 bw; 


nS >= ew; 


} 


B:n+2 ram := ¢; 


/*®* store overflow * / 


/* store sign bit */ 


/* 


td 


/* page 1 of YSOBEL.GAL */ 


Function: 
ysobel takes the sobel transform of the input image in the north - 
south direction (y direction). In other words, this function finds 
gradient of the image about a pixel in the north - south direction 
by implementing the algorithm: | 


f 


Ysobel = (D + 2E + F) = (J + 2K +L) 
where the 3x3 neighbdérhood around the pixel "H" is defined as: - 


DE F 
G HI 

JI K OL 

Parameters: 

Gray scale input image (i.e. unsigned 
magnitude number ) 

image B = Sobel transform of input image for the 


north - south direction (a 2's complement 
format number). 


image A 


Additional Memory Used: 


Execution Time: 


4-bit image 60 cycles 

8-bit image 104 cycles 

16-bit image 192 cycles 
Notes: _ 


1) Image B is three bits larger than image A. 


ysobel (A, B) 


image A, /* input image */ 
B: /* ouput image (y sobel of input) */ 
{ 
int i 


n = size(A); 

/* Compute first partial sum */ 
A:0 ew := ram c := O; 

A:0 ns := Fam ew := @e; 

B:O0 ram °:= smc := cy; 


for (i = 1; i < n; i++) 


{ 
A:i ew := ram; 
A:ins := ram ew := 6; 
B:i ram := sm Cc 3:2 CY; 
} 


Bin ram := sm c ;:= 0; /* store overflow bit */ 


“(4 


/* page 2 of YSOBEL.GAL */ 


/* Compute second partial sum */ 


for (1 = 0; i<n+til1; itt) 


{ 
B:i ew := ram; 
B:ins := ram ew := w; 
B:i ram := sm Cc i= oy? 
} 
B:n+l ram := sm c := Q; /* store overflow bit */ 


/* Compute difference */ 


for (i = 0; 1 < n+ 23 it+) 
{ 
B:i ns := ram; 
ns := s7 
B:ins :# ram ew := ns; 
ns := n; 
B:i ram := sm c := bw; 


} 
B:n+2 ram := ¢c} /*®* store sign bit */ 


